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Abstract 

We consider the dependence of collective flow on the nuclear surface thickness in a 
Boltzmann-Uehling-Uhlenbeck transport model of heavy ion collisions. Well defined 
surfaces are introduced by giving test particles a Gaussian density profile of constant 
width. Zeros of the flow excitation function are as much influenced by the surface 
thickness as the nuclear equation of state, and the dependence of this effect is 
understood in terms of a simple potential scattering model. Realistic calculations 
must also take into account medium effects for the nucleon-nucleon cross section, 
and impact parameter averaging. We find that balance energy scales with the mass 
number as A~y, where y has a numerical value between 0.35 and 0.5, depending on 
the assumptions about the in-medium nucleon-nucleon cross section. 



PACS Indices: 25.70.-z, 02.70.Lq, 21.65.+f 



Broadly speaking, nuclear collective flow in a heavy ion collision is the deflection 
of nuclear matter perpendicular to the beam axis during the course of the reaction. 
Experimentally, one observes that flow disappears at a well deflned beam energy 
P, the so called balance energy Esai, whose value depends on the system and 
impact parameter range being considered. These zeros in the flow excitation function 
were predicted by the Boltzmann-Uehling-Uhlenbeck (BUU) transport model @, 0], 
and an analysis of scale invariant quantities and may be understood as an 
overall cancellation of the attractive part of the mean fleld interaction with repulsive 
contributions from the mean fleld and coUisional kinetic pressure. Thus, as has been 
shown explicitly for BUU simulations in Ref. [Q], Esai is expected to depend on both 
the nuclear equation of state and the magnitude of the in-medium nucleon-nucleon 
cross section. By making a systematic study of the balance energy as a function of 
the nuclear mass, one therefore hopes to gain insight into these properties. However, 
other parameters might well influence the balance energy. In particular, we wish 
to investigate in this note the effect of flnite nuclear surface thicknesses and impact 
parameter variations on Esai- We begin by deflning the flow variable to be used 
here, and point out the importance of obtaining well-deflned nuclear surfaces that 
are independent of the grid size used to compute density gradients. We then show 
how the strong surface dependence of flow in a Vlasov simulation may be understood 
in terms of a simple potential scattering model. This dependence persists for full 
BUU calculations that include a non-zero collision integral. Lastly, we consider 
the effect on the balance energy when the mass number and impact parameter are 
varied. 

To analyze flow quantitatively in experiments one of the main problems that 
has to be addressed is the determination of the reaction plane (see, for example, 
Refs. In a model calculations, on the other hand, knowledge of the reaction 

plane immediately allows one to deflne a flow variable such as the average in-plane 
transverse momentum 
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where the weight = 1 or — 1 if the test particle i is emitted into the forward 
or backward center-of-mass hemispheres, respectively, and Pf is the transverse 
momentum of the test particle in the reaction plane. In this note, we shall refer to 
"flow" in the sense of the above equation. 

We begin by examining the way the local particle density p is calculated for the 
mean field dynamics in a Vlasov simulation. While the mean field is momentum 
dependent |T^-|T^, we shall consider here, for illustrative purposes, a Skyrme-like 



parametrization that is a function of the density alone: 

U{p)=A^ + B(P-y , (2) 
Po \PoJ 

and take values of the parameters A = —124 MeV, B = 70 MeV, and a = 2. This 
choice reproduces known nuclear matter properties, with a rather stiff compression 
modulus at saturation of Kq = 380 MeV. 

We choose to represent the nucleon phase space distribution function fifj-p^t) 
by an ensemble of test particles. If fif^p^t) is to satisfy the Vlasov equation, the 
equations of motion of a test particle with coordinates (r'j, Pi) are given by Hamilton's 
equations of motion with potential @): 

Pi = -^rU{p{fi)) and fi = ^!= = , (3) 



'mj, + Pi 

where m^v is the free nucleon mass. The local particle density p is often calculated on 
a grid, and the gradient obtained as a finite difference. In this procedure, each test 
particle counts a certain fraction towards the density of the cell it occupies, while 
neighboring cells receive a smaller contribution ||15[. This stabilizes the numerics, 
but also introduces a nuclear surface whose thickness is roughly given by the grid 
size As. This grid size must be larger than the maximum distance traversed by a 
test particle in one time step, but small enough to be able to compute the gradients 
in Eq. to sufficient accuracy. 

Instead of "smearing" a test particle in steps over only the nearest neighboring 
cells, one may choose to supply each test particle with, say, a Gaussian density 
profile of constant width. This also introduces a finite surface thickness, but one 



that is well-defined and independent of the grid size. Thus one may study the effect 
of varying the nuclear surface thickness without the external numerical constraints 
imposed on the choice of Ax. 

As an illustrative example, we consider flow in pure Vlasov dynamics as a 
function of the size of the grid for ^^^La on ^^^La collisions at a beam energy of 
800 MeV/nucleon and impact parameter b = 2.7 fm. Fig. 1 shows that (wPx) de- 
pends very strongly on the grid size, i.e., the nuclear surface (dotted line). The 
extrapolated value to zero grid size (zero surface thickness) is twice as large as the 
value at Ax = 1 fm, a frequently used grid size. This strong dependence per- 
sists at beam energies of 200 MeV/nucleon. At this energy, calculations with finite 
Ax > 0.5 fm gave an overall attraction, whereas the extrapolation to Ax = 
predicted positive flow. 

In Fig. 1, we also show the flow obtained with Gaussian test particle density 
profiles. The solid line represents the variation with Ax for a surface thickness of 
2 fm, while the dashed line is for a surface of 1 fm. Clearly, as long as the grid 
size is smaller than the Gaussian smearing width (i.e. the nuclear surface), the flow 
is independent of Ax. Also, the dotted line crosses the solid and dashed lines at 
Ax ~ 2 and 1 fm, respectively. This indicates that at least most of its rise with 
decreasing grid size is directly attributable to the changing surface thickness, and 
is not a "numerical artifact." We conclude that (wPx) depends rather strongly on 
the surface thickness, a thinner surface producing more than a thicker one. Of 
course, quantitatively the results obtained here are not reliable, since we have, for 
instance, ignored the momentum dependence and hard collisions, but they do show 
the important influence of the nuclear surface. 

The surface dependence of flow can be understood in terms of a simple potential 
scattering model |]14|. Assuming that the nuclei pass through each other without 
changing their shape in phase space, the centers of mass of the nuclei move according 
to the Hamilton function 
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where the potential V is given by 

V{R) = I V{p,ir))dr' , (5) 



M is the nuclear mass, Pi and P2 the center-of-mass momenta, and R the separation 
of the nuclei. In (^), V^pn) is the potential energy density corresponding to Eq. (|^) 

V^^^^^El^^P^, (6) 

To describe nuclei with a surface, we choose the density profile to be 



pR[r) 



1 + exp {\f-R/2\ -Rq) /a + 1 + exp (|f + R/2\ - i?o) /a 
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where Rq is the nuclear radius, and 4a the surface thickness. 

The potential (|^) is shown in Fig. 2 for various values of the parameter a. 
As expected, for an increasing surface thickness, the potential V decreases, and 
we expect a larger surface to produce less flow. This can be shown explicitly in 
the time evolution of the (wPx) (see Fig. 3). We find reasonable agreement of the 
potential scattering model with the test particle Vlasov calculation at both values of 
the surface thickness. Of course, differences are seen in the details, and are expected 
because of the crude assumptions made in Eq. (Q). For example, a Vlasov calculation 
shows that Esai occurs between 200 MeV and 300 MeV, depending on the surface, 
while the scattering model shows no zeros in the flow excitation function. 

We now focus on the disappearance of flow in a more realistic calculation, i.e., 
a full BUU simulation that includes collisions. For the reaction ^'^^La on ^^^La [|l| 
we find that the balance energy is shifted by ~ 10 MeV when the surface thickness 
increases from 1 fm to 2 fm. This is comparable to the shift expected when one 
changes from a stiff to a soft equation of state. For example, in ^°Ar on V reactions, 
Esai changes by only 8 MeV if the incompressibility is increased from Kq = 200 MeV 
to 380 MeV by adjusting the parameters in Eq. (|^) p. 

This becomes even more apparent in smaller systems such as ^^C+^^C, for which 
values of Esai were recently measured at the National Superconducting Cyclotron 
Laboratory |I6[ . Fig. 4 shows the flow obtained in simulations as a function of beam 



energy, for surface thicknesses of 1 fm (left panel) and 2 fm (right panel). The stiff 
EOS gives EBai ~ 140 MeV and 190 MeV, respectively, while the soft EOS yields 



^Bai ~ 170 MeV and 220 MeV, respectively. [|T^. The relative importance of the 
surface thickness is enhanced in smaller systems. Of course, the use of momentum 
dependent mean fields, a local Thomas-Fermi momentum space initialization, and 
values of the in-medium NN cross section will also influence the balance energy. 
For instance, it has been shown [Q] that BUU calculations with reduced NN cross 
sections result in higher values of the balance energy. 

Fig. 5 shows the balance energy as a function of the combined mass of the 
system. Experimental data (squares and diamond) suggests that the dependence of 
the balance energy on the combined mass. A, of the colliding nuclei follows a power 
law, 

EBai = xA-y , (8) 

where the exponent y has a numerical value of 0.33 ± 0.04. This power law depen- 
dence is reproduced by the BUU simulation (triangles and circles correspond to a 
soft and stiff EOS respectively) [|l^. In this calculation, we have parametrized the 
the nucleon-nucleon cross sections in terms of a least squares fit to the experimen- 
tal data of Ref. [O. The resulting values are somewhat different from the iso-spin 
averaged expressions used at higher beam energies. In particular, Cpp ^ (Tp„, and 
the cross sections are larger than the ones described in Ref. |]T5[. This results in a 
shift of the balance energy to lower values. 

In Fig. 6, we investigate the dependence of the balance energy on the value of 
the in-medium nucleon-nucleon cross section, where we look for medium corrections 
beyond the effect of the Pauli principle on the outgoing scattering states. In previous 
studies the free nucleon-nucleon cross section was multiplied with an overall constant 
scaling factor |^, ^ . However, this approach fails when one has collisions in low- 
density nuclear matter, where the in-medium cross section should approach its free- 
space value. A more realistic approach uses a Taylor expansion of the in-medium 
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cross section in the density variable: 



o-nn(v^,p) = (TNN(ys, 0) + p 



dp 



p=0 



= |^l + «i-^jaNN(yi,0) + ... (9) 

where we have introduced the dimensionless parameter ai, given by: 

d 

"1 = Po ^ (lnaNN(ys,p)} (10) 
op ^ J P=0 

In principle, ai is dependent on ^/s, but we have here - as a first approximation 
- taken ai as an energy- independent constant. Fig. 6 shows the mass dependence 
of the calculated balance energy as a function of different values of ai, where we 
have used a soft nuclear equation of state. It is clear that we obtain the best overall 
agreement with the experimental data for a value of ai = —0.2, correponding to a 
20% reduction of the nucleon-nucleon cross section at p = po. However, the power- 
law exponent y from Eq. (j^) depends on the choice of ai, with ai = (no medium 
modification) yielding the best agreement {y = 0.38 ± 0.05) with the experimental 
value of y = 0.32. For ai = —0.1 one extracts y = 0.42 ± 0.03, and for ai = —0.2 
we obtain y = 0.47 ± 0.03. That y increases with |q;i| is due to the fact that at 
higher beam energies higher densities are reached, and therfore the reduction of 
the in-medium cross section is stronger for the lighter systems, which have higher 
balance energies. This effect may, however, be at least partially compensated once 
one incorporates a more realistic energy dependence of ai. 

In addition to effects mentioned previously, the proportionality constant x (and, 
in principle, also the exponent y) in Eq. (P) also depends on the impact parameter 
b as can been seen in Fig. 7. For peripheral collisions, the contribution of the 
nucleon-nucleon collisions becomes relatively less important for the flow produc- 
tion than the mean fleld, because the repulsion generated by the nucleon-nucleon 
collisions is proportional to the overlap volume. We note that the calculations in 
Figs. 5 and 6 use only one impact parameter, obtained from the average value in 
the experiment, and proportionally scaled b for calculations at mass numbers A for 
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which no data exists. Of course, for a quantitative comparison with experiment a 
impact parameter weighted averaged flow has to be considered, but calculating Esai 
for several A and b is computationally prohibitive. 

In conclusion, we have shown that in Vlasov and BUU simulations reliable results 
for flow must take into account the finite thickness of the nuclear surface. This is 
best done by giving the test particles a Gaussian density profile, with a width that 
is larger than the grid size used to obtain the density gradients. The value of the 
nuclear surface thickness has a strong effect on the balance energy. Therefore, more 
quantitative BUU predictions have to not only take into account the equation of 
state and medium effects on the nucleon-nucleon cross section, but also proper initial 
conditions in phase space and impact parameter averaging. In addition, we find that 
a realistic variation of the in-medium nucleon-nucleon cross-section with density has 
a clear effect on the mass dependence of the balance energy. We find best overall 
agreement with the experimental data for cti = —0.2, where ai is defined in Eq. 
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Figure Captions 

Fig.l Flow (as defined in Eq. (JlD) versus grid size for pure Vlasov dynamics witli 
smearing over neigliboring cells only (dotted line). The dashed and solid lines show 
results if each test particle is given a Gaussian density profile of constant width. 
The dashed line corresponds to a nuclear surface thickness of 1 fm, the solid line to 
2 fm. All curves shown are for ^'^^La + La, E = 800 MeV/nucleon, and b = 2.7 fm. 

Fig. 2 Potential Eq. (^ as a function of the separation of the two nuclei. The 
radius Rq = 5.8 fm corresponds to ^'^^La. 

Fig. 3 The time evolution of flow for pure Vlasov dynamics (solid line), and the 
potential scattering model, Eq. (^) (dashes). The surface thicknesses are 1 fm (upper 
two curves) and 2 fm (lower two curves), and the system is as for Fig. 1. 

Fig. 4 Flow versus beam energy in the vicinity of the zeros of the excitation 
function. The system is ^^C C at an impact parameter 6 = 1.4 fm. The left 
panel is for a surface thickness of 1 fm, the right for 2 fm. Both panels show results 
for stiff (diamonds) and soft equations of state (plus signs). 

Fig. 5 The calculated values of Esai as a function of the mass of the system. Only 
symmetric systems are considered. Diamonds and circles correspond to a soft and 
stiff EOS, respectively. For comparison, experimental data from Ref. (squares) 
[^] are shown. The solid and dashed curves are power-law fits to the calculations 
and data. 

Fig. 6 Same as Fig. 5, but varying the in-medium cross section according to Eq. 
(^ with «! defined in Eq. (0). For all calculations, a soft equation of state was 
used. The lines represent power-law fits to the calculations and data. 

Fig. 7 The balance energy for CI + CI as a function of impact parameter. 
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